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Abstract 

We study the implications of minimal non-universal boundary conditions in the 
sfermion soft SUSY breaking (SSB) masses of mSUGRA. We impose asymptotic 
b-T Yukawa coupling unification and we resort to a parameterization of the devi- 
ation from universality in the SSB motivated by the multiplet structure of SU{5) 
GUT. A set of cosmo-phenomenological constraints, including the recent results 
from WMAP, determines the allowed parameter space of the models under consid- 
eration. We highlight a "new coannihilation corridor" where x~ bi and X ~ '^i ~ 
coannihilations significantly contribute to the reduction of the neutralino relic den- 
sity. 
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1 Introduction 



The minimal supersymmetric extension of the standard model (MSSM) and grand uni- 
fication are often regarded as the main ingredients of the physics beyond the standard 
model miSISl- The phenomenological implications of the MSSM and of grand unification 
have been investigated since decades, and, though direct evidence of such theories is still 
missing, a large set of increasingly stringent constraints has put important bounds on the 
parameter space of these so-called SUSY-GUTs. 

Nonetheless, the theoretical frameworks of SUSY-GUTs are countless, and typically char- 
acterized by a huge number of parameters. The common phenomenological practice is to 
make a certain number of (hopefully) theoretically motivated assumptions in order to deal 
with a reduced set of parameters. The first assumption one has to make in the context 
of SUSY is how to parameterize the mechanism of SUSY breaking. 

Here, we resort to one of the most widely studied such contexts, the one of the so called 
supergravity models. In this framework, SUSY is broken in a hidden sector, whose fields 
couple only gravitationally to the MSSM fields Under some further assumptions 
(gaugino and scalar universality), the soft SUSY breaking part of the Lagrangian of su- 
pergravity models is determined by four continuous parameters plus one sign (mSUGRA). 
In this paper, based on the results presented in we study a {minimal) deviation from 
universality in the scalar fermions sector of the theory, inspired by the multiplet structure 
of the simplest GUT, SU{5) (see sec l2.2|l . We impose to the model the constraint of b-r 
Yukawa coupling unification, which is a common prediction of a wide set of theories of 
grand unification, among which SU{5) and SO{10) GUTs. 

One of the virtues of i?-parity conserving SUSY models is to produce ideal candidates 
for cold dark matter jBJ , in the present case the lightest neutralino (a linear superposition 
of the neutral gauge and Higgs boson superpartners). The neutralino is, in fact, a color 
and electromagnetically neutral, weakly interacting massive particle (WIMP), as required 
to be a good cold dark matter candidate. The recent results from the WMAP satellite 
[7j, combined with other astrophysical data, determined with unprecedented accuracy the 
dark matter content range of the Universe. It has been shown [HI IH] that the cosmologi- 
cally allowed parameter space of mSUGRA is dramatically restricted by the considerably 
lowered upper bound on the neutralino relic density. Moreover, as recently pointed out in 
ref. j^, in cosmological scenarios involving quintessential fields, the relic density could 
undergo a further significant enhancement, even of several orders of magnitude. The 
necessity of efficient relic density suppression mechanisms is therefore by now an uncon- 
troversial point. 

One of the known mechanisms of neutralino relic density suppression occurs when the 
next-to-lightest supersymmetric particle (NLSP), is close in mass to the neutralino. In 
this case, the neutralino relic density is not only suppressed by neutralino-neutralino 
annihilations, but also by coannihilations with the NLSP, and indirectly also by the an- 
nihilations of the NLSP itself [HI 112 CHI- 

This mechanism can involve in mSUGRA various sparticles jHj, such as the lightest stau 
[TH| ITHt ITT] , stop [IHl UHl , or chargino [201 HE] • In the present paper we address the pos- 
sibility that minimal non-universality in the sfermion sector can produce "new" , unusual 
coannihilation partners. We find that this possibility is effectively realized, in what we 
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call a new coannihilation corridor, where the sbottom and the tau sneutrino play the role 
of NLSP. We show that this pattern is compatible with b-r exact ("top-down") Yukawa 
unification, and with all known phenomenological requirements. 

In the next section we introduce and motivate the model we analyse. We discuss the 
issue of b-T YU and of a GUT-motivated minimal sfermion SSB masses non-universality. 
We then show in secEJthe features of the resulting particle spectrum, which gives rise to 
"new coannihilation corridors" involving, besides the stau, the lightest sbottom and the 
tau sneutrino. secElis devoted to the description of the cosmo-phenomenological require- 
ments we apply. We demonstrate in particular that the /i > case is not compatible with 
b-r "top-down" YU. In secjHlwe describe the new coannihilation regions. After our final 
remarks, we give in the Appendix the complete list of the coannihilation processes involv- 
ing neutralino, sbottom, tau sneutrino and stau. Finally, an approximate treatment of the 
neutralino relic density contributions coming from sbottom-sbottom, sbottom-neutralino 
and stau-tau sneutrino coannihilations is provided. 

2 The Model 

2.1 b-r Yukawa unification 

One of the successful predictions of grand unified theories is the asymptotic unification 
of the third family Yukawa couplings Ij . The issue of Yukawa unification (YU) has been 
extensively studied, see e.g. |21t,22, In particular, in this paper we address the issue 
of YU of the bottom quark and of the tau lepton, which is a prediction of some of the 
minimal grand unification gauge group, such as SU{5). b-r YU is a consequence of the 
fact that the two particles belong to the same SU{5) multiplet, and therefore, at the scale 
of grand unification Mq^j, they are predicted to have the same Yukawa coupling. The 
experimental difference between rrir and rub is then mainly explained by two effects. First, 
the renormalization group (RG) running from Mqut to the electroweak scale drives the 
two masses to different values. Second, in the minimal supersymmetric standard model, 
the supersymmetric sparticles affect with different finite radiative corrections the values 
of the masses, in particular the one of the b quark 

Previous investigations of b-r YU include ref. |23 as regards non supersymmetric GUTs 
and ref. |2Sl EH EHl EH] , and more recently ref. [SHI EI] for the SUSY-GUT case. In 
particular, in jSU] the implications of the recent experimental, and theoretical, results on 
the muon anomalous magnetic moment and on the inclusive branching ratio b —>■ s'j were 
also taken into account, while in |j32j the neutralino relic density constrain was examined, 
in the context of gaugino non-universality. In ref. |33[ 13^ l85] the puzzle of neutrino 
masses and mixing has been tackled within the framework of b-r YU. 
A possible approach to b-r YU is of the "bottom- up" type [211 ED] • It consists in defining 
some parameter which evaluates the accuracy of YU, such as 

" hr{MavT) 

The procedure we take here is instead a "top-down" approach [SH EH EH] : for a given 
set of SUSY parameters we fix the value hT-^Mcuj.) = hb^Mqjjj.) requiring the resulting 
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rrir to be equal to its central experimental value. We then compute mf,(M^) through RG 
running and taking into account the SUSY corrections. A model giving a value of the 
6-quark mass lying outside the experimental range is ruled out. With this procedure, we 
perform exact h-r YU at the GUT scale, and directly check whether a given model can, 
or cannot, be compatible with it. 



2.2 Minimal sfermion non- universality 

The parameter space of the minimal supersymmetric extension of the standard model, 
in its most general form, includes more than a hundred parameters [SEIISZI- Therefore, 
it is commonly assumed that some underlying principle reduces the number of parame- 
ters appearing in the soft supersymmetry breaking (SSB) Lagrangian. In particular, in 
the case of gravity mediated supersymmetry breaking, one can theoretically motivate jlj 
the assumption that there exists, at some high energy scale Mx a common mass mg for 
all scalars as well as a common trilinear coupling term Aq for all SSB trilinear interac- 
tions. Moreover, in SUSY GUT scenarios, the additional assumption that the vacuum 
expectation value of the gauge kinetic function does not break the unifying gauge sym- 
metry yields a common mass M1/2 for all gauginos. One is then left with four parameters 
(mo, Aq, Ml/2, tan/?) and one sign (sign/i), which define the so-called constrained MSSM, 
or mSUGRA, parameter space. 

Much work has been done in the investigation of non-universality in the gaugino sector, 
see e.g. ref. jHElEnj- As regards the SSB scalar masses, it has been since long known 
that universality is not a consequence of the supergravity framework, but rather an addi- 
tional assumption jlU]. This justified an uprising interest in the possible consequences of 
non- universality in the scalar sector In particular, in an analysis of various 

possible deviations from universality in the SSB was carried out. 

In this paper we focus on a simple model exhibiting minimal non-universal sfermion masses 
(mNUSM) at the GUT scale. Our model is inspired by an SU{5) SUSY GUT where the 
scale of SSB universality Mx is higher than the GUT scale Mq^j, The RG evolution 
of the SSB from Mx down to Mq^-j, induces a pattern of non-universality in the sfermion 
sector dictated by the arrangement of the matter fields into the supermultiplets: 



L, D^) ^ 5 (1) 
Q, ^7^ i^) ^ 10 (2) 



This structure entails the following pattern of sfermion mass non-universality at the GUT 
scale: 

m\ = mjy = m| (3) 
^1 = ^11 = ^1 = ^lo (4) 

The running between Mx and Mqut will also produce two other effects. First, a typically 
large deviation from universality and a splitting between the up and down Higgs masses 
rriH-^ and rriH^ is generated at the GUT scale (a detailed study of non-universal Higgs 
masses is presented in jl2] and Second, a small splitting is also present between the 
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SSB masses of the two hghtest sfermion famihes and the third one. In the present paper we 
will however restrict to a phenomenological parameterization of sfermion non-universality, 
simply setting 
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We are therefore left with a single parameter which scans this "minimal", GUT- 
inspired deviation from universality in the sfermion sector. 

For our purposes, we let K vary between and 1: in this way we recover full universality 
(i.e. the CMSSM) for K = 1, while, for < 1, we can lower the spectrum of the 
sparticles belonging to the 5 multiplet. Hence, we generate a spectrum with significantly 
lower masses for the tau sneutrino and the lightest stau and sbottom. Whenever the 
masses of these sparticles are close to the neutralino mass, they can play an important 
role in coannihilation processes. 

Being inspired by a SUSY GUT SU{5) framework, it is natural to mention, within the 
proposed mNUSM model, the critical question of the proton decay j4S]. First, the non- 
universality pattern of mNUSM is not derived from a definite SU{5) GUT: it simply 
inherits from such theory a plausible asymptotic soft sfermion mass structure (w^io and 
m|) and the feature of b-r YU. Nonetheless, it has been shown that consistent SU{5) 
models exist [IH], where suitable structures for the leptoquarks Yukawa couplings Hqq, 
huE, and hqi drastically suppress the proton decay rate, even at large tan/5. The 
resulting proton life time is then well below the current experimental limits |37]. These 
consistent models are compatible with the present mNUSM in the soft SUSY breaking 
Lagrangian, though for computational ease we take here into account only the third 
generation Yukawa couplings ht and hf, = h^-. Hence, we conclude that mNUSM models, 
within an SU{b) framework, are viable and are not in contrast with the present constraints 
on the proton life time. 



3 Numerical procedure 

The mNUSM model we propose is defined by the following parameters: 

mo, Ao, Ml/2, tan/3, sign/i and i^. (8) 

We impose gauge and h-r Yukawa coupling unification at a GUT scale Mq^j, self-consistently 
determined by gauge coupling unification through two-loops SUSY renormalization group 
equations |1H], both for the gauge and for the Yukawa couplings, between Mqut and a 
common SUSY threshold Msusy — ^/m^^m^^ (ti^2 are the stop quark mass eigenvalues). 

At Ms 

usY we require radiative electroweak symmetry breaking, we evaluate the SUSY 
spectrum and calculate the SUSY corrections to the b and r masses j2H]- For the latter 
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we use the approximate formula of ref . [IH] ■ 
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m = mm[mi, m2j. 



(11) 



From MsusY to the running is continued via the SM one-loop RGBs. We use fixed 
values inO] for the running top quark mass mt{mt) = 166 GeV, for the running tau lepton 
mass mr{Mz) = 1.746 GeV and for as{Mz) = 0.1185, all fixed to their central experi- 
mental values. The asymptotic Yukawa couplings hT-^Ma^j.) = hb^MQ^j.) and hf^Ma^j.) 
are then consistently determined to get the correct top and tau masses, while the tree 
level m^f^'^'^ and the SUSY-corrected ml°" masses of the running bottom quark at Mz are 
outputs. 

The neutralino relic density is computed interfacing the output of the RGE running with 
the publicly available code micrOMEGAs [51j, which includes thermally averaged exact tree- 
level cross-sections of all possible (co-) annihilation processes, an appropriate treatment 
of poles and the one-loop QCD corrections to the Higgs coupling with the fermions. The 
output of micrOMEGAs also produces the relative contributions of any given final state to 
the reduction of the neutralino relic density. 

The direct and indirect detection rates are estimated through another publicly available 
numerical code, darkSUSY |52j . 

As regards the phenomenological constraints, the Higgs boson masses are calculated us- 
ing micrOMEGAs which incorporates the FeynHiggsFast [53^ code, where the SUSY 
contributions are calculated at two-loops. The inclusive BR{b — >■ S7) is again calculated 
with the current updated version of the micrOMEGAs code |33], where the SM contribu- 
tions are evaluated using the formalism of ref. j^SI and the charged Higgs boson SUSY 
contributions are computed including the next-to-leading order SUSY QCD re-summed 
corrections and the tan/3 enhanced contributions (see ref. ISEI). The SUSY contributions 
to the muon anomalous magnetic moment are directly calculated from the formulae 
of ref. |SZ1 and compared with the output of the micrOMEGAs code. 

4 The particle spectrum with mNUSM 

In minimal supergravity, especially after the very precise results from the WMAP satellite 
[7j, giving a considerably reduced upper limit on the CDM density of the Universe, the 
cosmologically allowed regions of the parameter space are strongly constrained [51in|. As 
pointed out in [5H1 EH] , naainly two mechanisms can suppress the neutralino relic density 
to sufficiently low values. The first one are coannihilations of the neutralino with the 
next-to-lightest sparticles (NLSP), which are effective whenever the mass of the latter 
lies within 10-20% of the neutralino mass |^ ITTl IT^ . The second mechanism is given 
by direct, rapid s-channel annihilation of the neutralino with the CP-odd Higgs boson A 
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[UT] , which takes place if ~ 2m^, i.e. if the channel is enhanced by pole effects^. In 
particular, we choose to focus here on the case of coannihilations, which in the present 
scenario of mNUSM are expected to exhibit a rich pattern. 

In order to understand the possible pattern of coannihilations emerging from the parame- 
terization of sfermion masses non-universality outlined in sec J2.2l we study the spectrum 
of the candidate NLSPs, namely the lightest sbottom and the tau sneutrino. The depen- 
dence of the respective masses of these two sparticles on the high energy SSB inputs can 
be parameterized as follows: 

a-,K'ml + {b-,ml + c~,Mlj^ + A^J (12) 
a^^K^ml + {hi,^ml + c^^M^^ + A^.^ (13) 

The contributions A^^ originate in part from the SU{2)l and U{1)y -D-term quartic 
interactions of the form (squark)^(Higgs)^ and (slepton)^(IIiggs)^, and are oc cos(2/?)M| 
|2j. The sbottom mass gets a further contribution in A^ arising from the LR off-diagonal 
elements of the mass matrix. The mixing terms are generated by the typically large 
values of A^, in its turn induced, even for Aq{Mg,jj.) = 0, by RG running, and by a 
tan /3-enhanced contribution oc ^rrih. In the case of the tau sneutrino, instead, a further, 
though small, contribution to Aj^^ analogously arises from Ar- 

We plot in figH] (a) and (b) the typical behaviors of the sparticle masses as functions 
of K, at fixed M1/2, mo and tan/3. Frame (a) shows the case of the mass of the sbot- 
tom, (we also plot the corresponding masses of the lightest neutralino and chargino). 
The high energy parameters are fixed at M1/2 = 1100 GeV, tan/5 = 38.0, Aq = and 
mo = 2850, 3300 and 3750 GeV. We clearly see from the figure that the behavior is, as 
expected, m^^ ^ ^/ a + f^K'^, where a < is the sum of the terms in parenthesis in eq. lfT^ 
and P = ai ml > 0. Increasing the value of mo induces higher negative values for a, as 
well as obviously higher values for /3. Notice in any case that for sufficiently low values of 
K the mass of the sbottom is driven below the mass of the neutralino and also to negative 
values. 

Frame (b) shows instead the mass of the tau sneutrino, at M1/2 = 1100 GeV, tan j3 = 38.0, 
Aq = and mo = 1350, 1650 and 1950 GeV. We see that also here rrij^^ ^ a/o + f3K'^, 
but in this case the interplay between bf^^ and Cj^^ generates either a positive (mo = 
1350, 1650 GeV) or negative (mo = 1950 GeV). The outcome is therefore that the 
sneutrino mass can be lowered towards the mass of the neutralino for low values of K, 
depending on mo; the typical range of mo for which this is possible is always lower than 
in the sbottom case. 

We carried out a thorough investigation of the possible coannihilation regions, and we 
found that a good parameter is represented by the ratio (mo/Mi/2) at fixed tan/5, M1/2 
and Aq. In fig El we plot the coannihilation corridors which we found scanning the plane 
(mo/Mi/2, fsT) at M1/2 = 1.1 TeV, tan/? = 38.0 and Aq = 0. Within the red solid lines 

^ Since here itia — mn, where H is the heaviest CP-even neutral Higgs boson, the condition ijia ~ 2m 
imphes pole effects also for H. However, since the CP quantum number of the exchanged Higgs boson 
must match that of the initial state, only A exchange contributes in the 5* wave, while H (and h) 
contribute to the P wave. This implies a suppression, in the thermal averaged cross section, of a factor 
3/a;_F « 0.1 0.2 (see the Appendix). 
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< 20%. In the lower part of the figure, which we indicate as 'Excluded Re- 

gion', the low value of mo implies that the stau becomes lighter than the neutralino, or 
even gets a negative (unphysical) mass. Very stringent bounds [02] indicate that the LSP 
has to be electrically and color neutral, and therefore this region is excluded. The strip 
above this excluded region represents a first coannihilation corridor, where the NLSP is 
the stau. We notice that this region survives up to = 1, i.e. fully universal boundary 
conditions. It actually represents a slice of the narrow band, in the (mo, M1/2) planes, 
which is cosmologically allowed thanks to neutralino-stau coannihilations (see e.g. fig.l 
and 2 of ref . 0). 

For values oi K < 0.5 we find a second, distinct, branch where the mass of the NLSP lies 
within 20% of the LSP mass. In the lower part of the branch (in fig|2lup to mo/Mi/2 ~ 3) 
the NLSP turns out to be the tau sneutrino, with the lightest stau which is quasi de- 
generate with it (say within few percent, see the discussion in sec J(i.2j) . Increasing mo 
and moving to the upper part of the branch, the lightest sbottom becomes the NLSP, 
while the tau sneutrino and the stau, still quasi degenerate, become by far heavier. In 
the upper left part of the figure, once again indicated as 'Excluded Region', either the 
sbottom or the tau sneutrino become lighter than the neutralino, or even get negative 
values for their masses (see figH}. In figOwe study the dependence of the shape of the 
coannihilation corridors on the various parameters which we fixed in figI21 In frame (a) we 
vary the values of M1/2, comparing the M1/2 = 1100 GeV case of fig|21with respectively 
Ml/2 = 1600 GeV and M1/2 = 600 GeV. We see that there is practically no significant 
dependence of the shape on the value of M1/2, and therefore we conclude that the chosen 
parameter mo/Mi/2 is good, being Mi/2-quasi independent. In frame (6) we vary instead 
the value of the universal trilinear coupling Aq. We plot again with a red solid line the 
Ao = case, as well as the Aq = 1 TeV and Aq = —1 TeV cases, always at tan/3 = 38.0 



7 



5 
4.5 

4 
3.5 

3 

OJ 

1 2.5 

O 

E 

2 
1.5 
1 

0.5 





0.2 0.4 0.6 0.8 1 

K 

Figure 2: Coannihilation regions in the {mQ/Mi/2i ^) P^<ine at — 1-1 TeV, tan/3 = 38.0 

and Aq = 0. Points within the red lines are characterized by — < 20%. In the top 

left and lower part of the figure the LSP is not a neutralino, while in the top right region no 
coannihilations take place between the neutralino and the sleptons. 



and Ml/2 = 1100 GeV. Once again we do not see any significant effect, apart from a 
common sliift of tlie lower coanniliilation brancfi upwards and of tlie upper downwards. 
Tliese sfiifts can be traced back to tlie effect of tlie off diagonal Aq (sign-independent) 
entries in the sfermion mass matrices. 

FiglU illustrates the dependence on tan/5, highlighting this time a significant effect on the 
upper branch: increasing tan/3 yields higher values for K, i.e. the branch is moved to the 
right, and vice- versa. This effect can be qualitatively understood from the approximate 
expression for the mass eigenvalues of the relevant sfermions, eq. fjl2|) . where an increase in 
tan (3 is compensated by an increase in the effective boundary value of the scalar masses, 
tuned by the parameter K. For instance, Aj,^ contains a negatively contributing term 
oc tan jSrribfi, which is compensated, when the value of mo and of is fixed, by an in- 
crease of K . 



5 Cosmo-phenomenological bounds 

In this paper we apply two classes of constraints: on the one hand the cosmological 
bounds coming from the limits on the cold dark matter content of the Universe and from 
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Figure 3: Coannihilation regions in the plane at tan/3 = 38.0 and Aq = 0. In 

frame (a) the dependence on M^j^ is studied at Aq = 0. In frame (h) two different values for 
the trilinear coupling = =b 1 TeV are plotted at = 1100 GeV. 
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Figure 4: Coannihilation regions in the {mQ/Mii2,K) plane at = 1-1 TeV and Aq = 0. 

The red solid lines indicate the coannihilation corridor for tan /3 = 38.0, while the green dashed 
and the blue dotted lines respectively tan /3 = 34, 42. 



direct and indirect neutralino searches; on the other hand we impose the most stringent 
"accelerator" constraints, such as the inclusive BR{b —* 57) and the Higgs boson mass. In 
this second class we also include the bound on the 6-quark mass, direct sparticle searches 
and a conservative approach jHS] to the constraint coming from the SUSY corrections to 
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the muon anomalous magnetic moment 5a^. 

For illustrative purposes, we show the behavior of three benchmark cases, pertaining the 
three regions of coannihilations highlighted in the previous section. Namely, we choose 
three representative values of K and require the mass splitting between the neutralino 
and the NLSP to be 10%. We then scan the parameter space, setting Aq = and fi < 
for simplicity and varying m^. The features of the three cases considered are summarized 
in Tabffl 



5.1 Neutralino relic density facing WMAP results 

Supersymmetric models with conserved R parity generate ideal candidates for cold dark 
matter, as, in the present case, the lightest neutralino. It is therefore natural to require, 
besides the fulfillment of the other phenomenological constraints, that the cosmological 
relic density of the neutralinos lies within the bounds indicated by cosmology. In partic- 
ular, the recent results from the WMAP satellite [7], combined via a global fit procedure 
with other astrophysical data (including other CMB experiments, LSS surveys and the 
Lya data), give a compelling bound on the cosmological relic density of Cold Dark Matter 



n /,2 - n 11 96+0-00805 

We take here the 2-a range, and we require that ^xh'^ < 0.1287. Strictly speaking, the 
lower bound cannot be directly imposed, if we suppose that the neutralinos are not the 
only contributors to the cold dark matter of the Universe. For illustrative purposes, in 
fig^^we will plot also the lower bound on the relic density. 

In figlHlwe show instead the behavior of the neutralino relic density as a function of rriy. for 
the three benchmark cases. We see that in the case of tau sneutrino the coannihilations 
only weakly contribute to the reduction of the relic density (see sec lfi.2j) . which, as ex- 
pected, quadratically diverges with m^. In the case of the sbottom coannihilation, SUSY 
QCD effectively enhance the relic density suppression, which nonetheless still exhibits a 
divergent behavior. In the case of the stau, instead, we notice how the interplay between 
coannihilation and direct rapid annihilation through the A-po\e s-channel can drastically 
reduce the relic density. In the dip located between 600 and 700 GeV, in fact, the split- 
ting between the mass of the CP-odd Higgs boson and twice the mass of the neutralino 

— ^ — \ — ^ ~ 3%, and therefore direct pole annihilations are extremely efficient, leading 
to viable values of ^^Jn? for rather high ^ 700 GeV. 



5.2 Direct and indirect WIMP searches 

As pointed out in fHS], non-universality in the sfermion sector should not give rise to 
substantial modifications to o^^^'^ In fact, owing to the small Yukawa couplings of 
the lightest generation of quarks, NUSM do not greatly affect the masses of the up and 
down squarks, leading to substantially unchanged values for the direct annihilation cross 

sections, sensitive to processes like X <? X 9- Further, we find that Higgs exchange 
processes do not give rise, in the mNUSM models, to sizeable contributions, as compared 
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Table 1: The three benchmark scenarios described in the text 
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Figure 5: Neutralino relic density for the three benchmark cases of Tab^ The upper and lower 
bounds on ^cDMh"^ are the 95% C.L. from WMAP global fit 



with the fully universal case. In figEl (ct) we scatter plot the spin-independent cross sec- 
tion CF^^^q as a function of for K = 1.0 (red +, full universality case), K = 0.35 (green 
x) and K = 0.1 (blue *). We randomly scan the parameter space for a given value of K 
with 50 < Ml/2 < 1000 GeV, 20 < mo < 3000 GeV and -500 < Aq < 500 GeV, fixing 
tan/3 = 38.0, and requiring the fulfillment of the phenomenological constraints. We also 
include the sensitivities of present and future direct detection experiments [MI. Notice 
that some mNUSM models lie within the reach of future direct detection experiments. As 
far as indirect detection (e.g. the muon flux from the Sun) is concerned, non-universality 
do not significantly affect the detection rates, though t-channel sfermions exchange is 
enhanced by lighter third generation sfermions for i^' < 1, as we can see in figEl {b): 
for smaller values of K we obtain larger muon fluxes from the Sun. Present and future 
indirect detection experiments are however sensitive to muon fluxes well above what 
we obtain for mNUSM models. 

To summarize, we find, both for direct and indirect detection, that the mNUSM sce- 
nario produces detection rates in the same range of the CMSSM [HE], therefore (minimal) 
sfermion non-universality can hardly be inferred from WIMP searches. Further, present 
experimental bounds do not constrain the models under consideration. 
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Figure 6: Direct (a) and indirect (h) detection rates for K = 1.0 (red +, full universality), 
K = 0.35 (green x) and K = 0.1 (blue *). The points refer to 50 < M1/2 < 1000 GeV, 
20 < mo < 3000 GeV and -500 < Aq < 500 GeV, at tan/3 = 38.0. In frame (a) we plot o-|:!^J^ 
and the sensitivity bounds of present and future direct detection experiments ^641- In frame (b) 
we report the neutralino annihilation induced muon flux from the Sun, in Km.~'^y~^ , and we 
compare again with current and future indirect detection experiments sensitivities 155^ . 



5.3 SUSY corrections to the 6-quark mass 

It has been already pointed out, see e.g. [HUl EI], that the requirement of b-r Yukawa 
unification favors negative^ values of fi. We recall that sign/i is one of the parameters 
included both in the CMSSM and in our mNUSM model. In this section we demonstrate 
that, even with mNUSM, /i > is not compatible with YU. We then discuss the fi < 
case which allows, in a suitable tan jS range, the fulfillment of b-r YU. 
The main problem of b-r YU with > is that one typically obtains a tree level mass for 
the b quark which is close to the experimental upper bound, and has to add on top of it 
large positive SUSY corrections (eq. ljlfij) ). which drive ml°" outside the experimental range 
(or, the other way round, hh far away from hr, in the bottom-up approach). We impose 
b-T YU at the GUT scale, we fix hT-{MauT) = hb{Mc^j.) from the properly corrected and 
RG evolved m^(M^), obtaining as outputs the tree level ■m]f^'^ and the SUSY-corrected 
^corr jYiasses of the running bottom quark at Mz- We compare these numbers with the 
appropriately evolved 6-quark pole mass up to the Mz scale, with as{Mz) — 0.1185, 
following the procedure of ref . IHS ED] : 



nibimb) = 4.25 ± 0.3 GeV mb{Mz) 
^We use here the standard sign conventions of ref. [HZI 
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= 2.88 ±0.2 GeV 



(15) 



The largest SUSY corrections arise from sbottom-gluino and stop-chargino loops, frozen 
at the MsusY scale |21l IIHI EI] • They are non- decoupling effects because one gets a finite 
contribution even in the infinite sparticle mass limit, and they can be cast in the following 
approximate form: 

^ ^ ^M3/i/K,m|,M|) tan/3, (16) 
^"^^ ~ ^* /(m?,m|,/x2) tan/3, (17) 



I{x,y,z) = 



167r2^ ' ^ 

xy \n{x/y) + xz ln(z/ x) + yz ln(y/ z) 
{x - y){y - z){x - z) 



Unless the trilinear coupling At is very large, the gluino loop typically dominate (an 
exception is investigated in ref. [72]) and the sign of the SUSY contribution is given by 
the sign of M^fi. Therefore, since we assume here gaugino universality, this implies that 
b-T YU is favored in the fi < case. 

To numerically quantify this statement, we study the behavior of m^f^'^'^ and m^°" varying 
the parameters of the mNUSM model. In fig[7j (a) we fix Aq = 0, tan/3 = 38.0, mo = 
1000 GeV and M1/2 = 1100 GeV, and analyse the dependence on K. Wee see that 
the tree level value ml^'^^ is roughly constant, while the SUSY corrections decrease as K 
increases. Both remain however well above the experimental upper bound. This can be 
understood from eq. (jlfij) . since increasing K means increasing m^^, with a fixed value for 
M3 and, roughly, for /i and mi^, and the function I{x,y,z), at fixed y and z is inversely 
proportional to x = . Therefore, subsequently, we concentrate on the universal K = 1 
case, in order to check whether a parameter space allowing for "top-down" b-r YU exists 
or not. 

The second step is to study the dependence of the 6-quark mass on the parameter mo (figlZI 
(b)). We take here Aq = 0, tan/3 = 38.0, M1/2 = 1100 GeV and K = 1, and we notice, as 
expected, that the size of the corrections decreases with increasing mg. This is explained 
on the one hand by the fact that from the radiative EWSB condition the value of /i^ is 
slightly decreased by the increase of mo |2] , and on the other hand because increasing the 
value of mjj^ ^ leads to a decrease of the function I. In figEI (c) we take instead Aq = 0, 
tan/3 = 38.0, mo = 2000 GeV and K = 1, and vary Mi/2- As can be easily inferred 
from eq. ()16|) . increasing M1/2 leads to an increase both in M3 and in /x. In conclusion, 
the candidate parameter space for b-r YU for /i > is at low M1/2 and high mo values. 
We choose therefore two trial values, M1/2 = 300 GeV, mo = 2000 GeV and we show 
our results in figlHl As readily seen from eq. (fTH|) . we find that the SUSY contributions 
grow with tan j3. We notice however that the tree level mass strongly decreases with tan /3, 
owing to the fact that the positive SUSY contributions to rrir (see eq. (jHl)) imply a smaller 
value for the asymptotic common b-r Yukawa coupling. The overall conspiracy of these 
two effects is to maintain the corrected b quark mass well above the experimental upper 
bound. This conclusion is further confirmed investigating extreme values of (M1/2, "^o) 
for high tan/3, always in the universal K = 1 case, and resorting also to nonzero values 
of Ao. 
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Figure 7: Tree level and SUSY corrected values of the b-quark mass at tan /? = 38.0 and Aq = 0. 
In frame (a) the parameter K is varied at fixed M1/2 = 1100 GeV and rriQ = 1000 GeV. In 
frame (b) tuq is varied at fixed M1/2 = 1100 GeV and K = 1. Frame (c) shows the dependence 
on Ml/2 at niQ = 2000 GeV and K = 1. 



To sum up, we demonstrated that, due to the SUSY corrections to the 6-quark mass, top- 
down b-T YU is excluded, both with universal and with minimal non-universal sfermion 
masses, in the case /i > 0. 



5.3.1 The < case 

In the fi < case the SUSY contributions to the 6-quark mass are negative, and therefore 
conspire to bring the tree level mass dictated by b-r YU within the experimental range 
ffT^ . We plot in figEl (a) the tree level and SUSY-corrected values of the 6-quark mass 
for the three benchmark cases of TabO We notice that the size of the SUSY corrections 
decrease when m^. is increased, because of the mentioned behavior of the function I{x, y, z) 
of eq. p(j|l . which is inverse-proportional to its arguments. In the case of the sbottom 
coannihilation (green dashed line) the effect is enhanced by the large size of uiq which 
directly reduces, through the suppression of /x, due to the EWSB condition, the size of 
the SUSY contributions. 



5.4 The inclusive branching ratio 6 ^ 57 

We construct the 95% CL. range on the BR{b 57) starting from the recent experimen- 
tal data of ref. and properly combining the experimental and theoretical uncertainties. 
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Figure 8: Tree level and SUSY corrected values of the b-quark mass at M1/2 = 300 GeV, mo 
2000 GeV, Ao = and K = I. 
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Figure 9: Tree level and SUSY corrected b-quark masses (a) and the BR{b — > 57) for the three 
benchmark scenarios of Tab^\ The upper and lower bounds are those of eq. U5\) for frame (a) 
and those of eq. il9\) for frame (b). 



The resulting bound is 



1.9 X 10"^ < BR{b 57) < 4.6 x 10"^ (19) 
15 



We calculate the BR{b — > sj) using the current updated version of micrOMEGAs [54 . The 
SM contribution is calculated through the formulae of ref. while the SUSY corrections 
through the ones of ref. [SEl- We can estimate the dependence of the SUSY corrections 
from the approximate formulae of ref. [731 EH IZS|- They are given, respectively for the 
charged Higgs and for the chargino exchanges, by 
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We notice that the relative sign of the two contributions of eq. ()2()j) and ()2H1 is given 
by sign(At/i) = — sign(M3yu), (the latter equality being valid for large top and bottom 
Yukawa couplings, as in the present case). Moreover, the SM contribution has the same 
sign as the charged Higgs one. All the contributions decrease as m<^ increases, therefore 
we expect to draw a lower bound on from the lower bound on BR{b s-y) in 
the case M^fi > and from the upper bound of eq. ()19|1 for M^fi < 0, which is the 
present case. We notice that the lower bound on becomes more restrictive as tan jS is 
increased, as can be read out from eq. (PT|) . 

In figini (b) we plot the results we get for the three benchmark cases of Tabffl We can see 
that as the SUSY spectrum becomes heavier (we recall that if the sbottom is the NLSP 
mo/Mi/2 ^ 3, see figEI) the lower bound on weakens, as can be understood from eq. 
()2()|) and ()2H1 and from the explicit form of f^^ and . 



5.5 The Higgs bosons masses 

We take in our analysis the 95% CL. LEP bound [70] on the lightest CP-even neutral 
Higgs boson mass 

TUh > 114.3 GeV, (24) 

which gives a lower bound on the m^. Two loop corrections are taken into account, as 
described in secj^l In figUHl (o) we plot the results for the Higgs mass we obtain for the 
three benchmark scenarios. 



5.6 Direct accelerators sparticle searches 

We impose the current direct accelerator search limits on the sparticle masses respec- 
tively > 104 GeV, > 100 GeV for / = ii,biJ^,i), rUg > 300 GeV, Trig,, > 260 

GeV for q = u,d, s, c. We always find the constraints from BR{b — > S'-f) and from rrih 
more restrictive than the direct sparticle searches. Therefore we do not plot these bounds 
in the panels of sec J6.1l and 16.21 
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Figure 10: rrih (a) and 6afi^^'^ (b) for the three benchmark scenarios of Tab^\ The lower bound 
of frame (a) is given in eq. \24]j , while in frame (b) the 2-ct and 5-0" bounds are taken from eg. 



5.7 The muon anomalous magnetic moment 

The deviation 5a^ of the muon anomalous magnetic moment from its predicted value 
in the Standard Model can be interpreted as arising from SUSY contributions, 6a^^^^ , 
mainly given by neutralino-smuon (a) and chargino-sneutrino (6) loops. The quantity 
5a^^^^, which we numerically compute with the micrOMEGAs package, can be estimated 
through the formulae of ref . , assuming a common value for the SUSY sparticles masses 

'^SUSY 

6ar- ^ sign{M,^.)^^-^{59l + 9!). (25) 

iyz,/i iii'susY 

On the experimental side, the BNL E821 experiment recently delivered a high precision 
measurement (0.7 ppm) of a^^^ = 11659203(8) x 10~^°. The theoretical computation of the 
SM prediction is plagued by the problem of estimating the hadronic vacuum polarization 
contribution. In particular, there is a persisting discrepancy between the calculations 
based on the r decay data and those based on low-energy e^e" data. Recent evaluations 
|77| 1751 17^ give the following range for the deviation of the SM value of from the 
experimental one: 

^exp_^sM _ g ^ ^^ 7) X 10^1° (e+ e") (26) 

^exv _ ^SM ^ ^^Q_4 ^ ^ ^Q-IO ^gp^y) (27) 

Following the hues of ref [HSl we decided to take a conservative approach to the problem 
of choosing which bound should be culled from the tantalizing results of eq. 1^ . 
In what follows we will therefore just indicate the region determined by the 2-a and 5-cr 
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range for the e"'"e~ based approach and for the r decay based approach, but we will not 
use this constraint to derive (lower) bounds on m^. 

In figEl (&), where we present the results we get for the three benchmark scenarios, we 
indicate only the r decay based limit, since 5a^^^^ , being < 0, is negative (see eq.(P^). 
In passing, we remark that it has been recently claimed that the result from e+e^ 
data is less clean than the one from r-decay (07|) . This would be due to intereferences 
of isoscalar / = mesons, not produced in r-decay, with vector mesons [ZO]- Therefore, 
theories requiring a negative sign of /x, as those with exact h — r YU, seem to be no longer 
disfavored by the 5a^ constraint. 

Formula (^3)1 allows us to interpret the results we show in fig^l(6) and in the figures of the 
next section. First, we see that as the SUSY spectrum becomes heavier (sbottom NLSP 
case) the corrections become smaller, as emerging from ()25|) . Second, we can see that 
^qSusy (jgcrggLses with increasing m^, for the same region as above, and that it increases 
with tan /3. 

6 The new coannihilation corridor 

In the yU < case the SUSY corrections to the 6-quark mass are negative, and therefore 
can naturally drive, for suitable values of tan/5, the corrected 7jff^"{Mz) within the ex- 
perimental range. 

In this section we focus on the extended sfermion coannihilatoin modes allowed by the 
particle spectrum of mNUSM. As pointed out in secHl the boundary conditions given by 
mNUSM produce a new "coannihilation branch", extending from K = Q xvp io K ^ 0.5 
and for values of mo greater than cMi/2, with c of o(l), increasing with decreasing val- 
ues of tan (3. In the upper part of the branch, typically for mo/Mi/2 ^ 3, the NLSP turns 
out to be the lightest sbottom, while for c < mo/Mi/2 < 3 we find ~ ~ m^^. The 
following two subsections are devoted to the analysis of the coannihilations processes x-^i 
and x-^T-T'ii taking into account the constraints described in the previous section. The 
final two subsections deal respectively with the determination of the tan 13 range allowed 
by h-T YU, and with the question of fine-tuning in mNUSM models. 

6.1 X~^i coannihilations 

Being a strong-interacting particle, we find that sbottom annihilation and coannihilations 
channels are very efficient, and substantially contribute to the neutralino relic density 
suppression. 

We plot in fig^2and^l(a) and (6) the cosmologically allowed regions for three represen- 
tative values of tan (3 = 34, 38 and 42 for a fixed value oi K = 0.35. The y axis represents 
the relative percent splitting between the neutralino mass and the next-to-lightest SUSY 
particle, in this case the lightest sbottom: 

A,,,, = (28) 

The relevant phenomenological constraints described in the previous sections determine 
lower, and sometimes upper, limits on m^. They result in almost vertical lines, since they 
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typically weakly depend on mo, which moreover very little varies in the plotted regions. 
They are instead fixed by M1/2 (and therefore, in the plots, by m^^). The regions allowed 
by a given constraint lie at the right of the respective lines. Points below the red lines are 
characterized by relic densities which fall within the 2-a range of eq. ()14|) . We emphasize 
that the hmits on and Ajvlsp we find here cannot be regarded as being representa- 
tive of the full parameter space: they show instead, for a few benchmark cases, how the 
interplay of the cosmo-phenomenological bounds described in secE] constrain the mass 
spectrum (here parameterized by rriy. and Ajvlsp) of the models under consideration. 
In fig^^we plot with a magenta dotted line also the lower bound on fi^DM^^- We see 
that the most stringent bound, for tan/3 < 38, is provided by the BR{b —>■ 57), while for 
higher values of tan f3 the 6-quark mass constraint becomes more restrictive. As far as the 
upper bound on is concerned, we find that high values of tan/3 allow to extend 
up to 1.5 TeV (figII21 (b)). For smaller tan/3 the upper bound on is instead fixed by 
the upper bound on the 6-quark mass: the SUSY corrections to (see eq.®) drive the 
common b-r Yukawa coupling to higher values, therefore generating a higher tree level 
6-quark mass which is not compensated by the (smaller) negative SUSY corrections to 
rrih (eq. (|T6|l ). The inverse mechanism pushes the lower bound on to higher values 
(^ 0.5 TeV) in the case of tan/3 = 42. We notice that the region determined by the 2-a 
5a^ range from the r decay approach always covers the great part (or the totality) of the 
allowed parameter space. 

In figEl we studied the sensitivity of our results to a non-zero value of the common 
trilinear coupling Aq at the GUT scale, for tan/3 = 38.0, A^^ ^ and K = 0.35. We 




Figure 11: Cosmologically allowed parameter space in the shottom coannihilation region, at 
K = 0.35 and tan/3 = 38.0. The region below the red solid line has VLyJi^ < 0.1287. The 
magenta line shows the putative lower bound on the neutralino relic density \14\ )- The blue 
dotted nearly vertical lines indicate the phenomenological constraints described in the text, while 
the green dashed line stands for the 2-a bound on the muon anomalous magnetic moment. 
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Figure 12: Cosmo-phenomenologically allowed parameter space in the sbottom Coannihilation 
region at K = 0.35, tan/3 = 34.0 (a) and tan/3 = 42.0 (h). 
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Figure 13: The determination of the lower hound on in the case of non zero Aq, for 
tan /3 = 38.0, A^^ ~ and K = 0.35. The allowed region is in the upper right part of the figure, 
as indicated. 



notice that the lower hmit on is fixed by the Aq = case. As for the cosmological 
bound on Qy-h"^, the variation of Aq leads to very little modifications of the upper bounds 
on TJiy. displayed in fig^2 and ^1 we therefore conclude that the only relevant change 
is the type of constraint that, depending on the sign and on the absolute value of Aq 
determines the lower bound on m^. 

As far as the coannihilation channels are concerned, the case of the sbottom is character- 
ized by a rather simple pattern, clearly dominated by strong interaction processes. We 
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find that for sbottom masses quasi-degenerate witli tfie neutralino mass, tfie neutralino 
pair anniliilation rate is very low (less than few percent). The dominant channels concern 
instead neutralino-sbottom coannihilations into gluon-6 quark (up to 10%) and sbottom- 
sbottom annihilations into a couple of b quarks (up to 15%) or into a couple of gluons 
(up to 80%). We give in the Appendix fsec ]A.2|) an approximate analytical treatment of 
these three most relevant processes. 

6.2 X — T^i — i^T coannihilations 

In the lower of part of the new coannihilation branch of fig for (mo/Mi/2) ^ 3 and 
< K < 0.25, the NLSP switches from the sbottom to the tau sneutrino. In this region 
the lightest stau is always heavier than the sneutrino, but the relative mass splitting is 
within few percent. This small splitting results from the combination of two RG effects: 

ml - ml ^ Al^iAr) - qcos(2/5)M|. (29) 

The first contribution in eq. ()29|l comes from the mixing between tl and f^j, which de- 
pends on Ar, while the second term stems from the mentioned D-teim quartic interac- 
tion, and the coefficient cj ~ 0.8 Even though = at the GUT scale, RG effects 
can drive At-^Ms^sy) to large values (for the case of figO we obtain typical values for 
At- ^ 0.5 TeV) and thus entail a non trivial LR mixing A£^. In the case of figOwe get 
ml — m? ^ 10 GeV. 

In fig^^we plot the cosmo-phenomenologically allowed region at K = 0.1 and tan/3 = 
38.0, in analogy with figHU We also show the "super-conservative" 5-a bound on 6a^, 
while the m^ bounds lie outside the depicted range. 

figim and El allow to compare the efficiency of coannihilation processes involving a 
strongly interacting sparticle (the sbottom) and a weakly, or electromagnetically interact- 
ing one (the tau sneutrino and the stau). We see that in the second case the maximum 
mass splitting between the NLSP and the neutralino is 2%, while in the first one, at the 
same m^, the cosmological bound allows a mass splitting up to ~ 8%. In figHSlwe plot 
the allowed regions at tan/3 = 34 (a) and 42 (6). We notice in all cases that slepton coan- 
nihilations constrain m^ to less than approx. 0.5 TeV, while in the case of the sbottom 
the bound is three times higher. 

Remarkably, we see in figEl and ^1 that the 2-a bound on (5a^ coming from r de- 
cay data is always fulfilled in the parameter space regions allowed by the other cosmo- 
phenomenological constraints. 

As regards non- zero values for Aq, we draw in the present case the same conclusions as 
in the preceding section (see fig lT!^ : the type of phenomenological constraint giving the 
lower bound on m^ in general changes, but the lowest value of m^ is determined by the 
Aq = case. 

The coannihilation pattern is in this case by far more complicated than in the sbottom 
case. We show in fig^Ja typical situation for the (percent) contributions of the possible 
coannihilating initial sparticles, and a detail of the most relevant final states, taken at 
m^ ~ m/>^ and tan/3 = 38. This pattern is however rather dependent on the tan/3 value 
and on the relative mass splitting. The case of stau-tau sneutrino coannihilations is fur- 
ther discussed in the Appendix. 
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6.3 The tan/3 range 



From eq. (|T^ we notice that the size of the corrections (hnearly) grows with tan (3, 
therefore we expect that the upper bound on mh{Mz) determine the lower bound on 
tan f] and vice-versa. As emerging from the figures of sec 16.11 and l(i.2| the constraint 
concurs with the other cosmo-phenomenological bounds, entails the determination of the 
allowed range of tan/3 fsec J6.3|) . We determine the lowest limit on tan/? combining the 
most restrictive phenomenological constraint, which turns out to be the BR{b — > 57) and 
which gives the lower bound on m™™ = m^, with the upper bound on ml°^^, which in 
its turn gives the upper bound m™^^. Requiring m™™ = m™^^ unambiguously gives the 
lowest allowed value of tan (3. In order to find the upper limit on tan j3, we notice that the 
bound on rrib is weaker than the constraint coming from BR{b —>■ sj) in determining the 
lowest value m™° (see e.g. figHH (&) and figlISl(&)). In fact, the BR{b — > S7) constraint 
becomes more and more stringent as one increases tan/?. On the other hand, m™^^ is 
this time fixed by the cosmological constraint on the maximum allowed neutralino relic 
density (see sec 15. II and figII21(&) andlTHl (&)). In figHJIwe depict the allowed tan/3 range 
in the (tan /3, A^vlsp) plane for three benchmark K values, respectively K = 0.8, 0.4 
and 0.1. Points below the lines are cosmologically and phenomenologically viable. As 
emerging from figEl in the case of stau NLSP the coannihilation region can be connected 
to the A pole region, hence the allowed Ajvlsp becomes large and no more meaningful 
to test the efficiency of coannihilation processes. FigEI also highlights that sbottom 
coannihilations allow a larger mass range for the coannihilating particle, since they involve 
strong interaction processes: the maximum Ajv^sp extends in this case up to ~ 13%, while 
in the case of x ~ t^i ~ i^r coannihilations its maximum value is around 3%. In the high 
tan P tail of figHZl as easily understood, x ~ ^1 coannihilations are slightly less efficient 
than X — fi — z/,-. Though our analysis was limited to benchmark K values, we can conclude 
from figEI that in the fi < case the allowed tan /? range is 

31 < tan/5 < 45, /i < 0. (30) 

We checked that the range given in ()30|1 holds also for Aq 0. We find in fact that in 
this case an analogous procedure of determination of (tan/5)™^^' ™™ yields weaker lim- 
its. What happens at 7^ is mainly that the phenomenological bounds tend to push 
m™™' to higher values, therefore (tan/3)™'^^ is lowered, while (tan/5)™™ is left substan- 
tially unchanged (see e.g. fig lT^ . 



6.4 Fine tuning 

The question of the quantification of fine-tuning in SUSY models for dark matter, be- 
sides involving the naturalness of the models, is also related to the possibility of reliably 
predicting the relic density ^x^"^ from accelerator measurements of the SUSY input pa- 
rameters jHOl • We emphasize that the present mNUSM models can hardly be regarded as 
realistic SUSY models, being conceived in order to address the issue of extended sfermion 
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Figure 14: Cosmologically allowed parameter space in the tau sneutrino-stau coannihilation 
region, at K = 0.1 and tan (3 = 38.0. The scale of the axis, as well as the notation, are the same 
as in fiallll 
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Figure 15: Cosmo-phenomenologically allowed parameter space in the tau sneutrino-stau Coan- 
nihilation region at K = 0.1, tan/? = 34.0 (a) and tan/? = 42.0 (h). 



coannihilation processes in the presence of non-universal sfermion masses boundary con- 
ditions. Nevertheless it is natural to face here the problem of the fine-tuning related to 
the extra parameter K appearing in mNUSM models. In fact, as emerging from fig^and 
fig 121 the new coannihilation regions appearing at large mo for K < 0.5 are expected to 
generate abrupt changes in the neutralino relic density for small variations of the mNUSM 



23 




XX XV X'f VV VT TT 



Figure 16: ^ typical pattern of the relative contribution of coannihilation processes in the tau 
sneutrino coannihilation region. The plot refers to the case tan/3 = 38, = and ~ rri;^^. 
The upper gray part of the columns represents other contributing channels. 
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1 K=0.8, NLSP=Stau 






1 K=0.4, NLSP=Sbottom 






1 K=0.1, NLSP=Tau Sneutrino 
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1 1 1 1 ''i"— j_ 



30 32 34 36 38 40 42 44 

Tan p 

Fi gure 17: The allowed (tan P, Aj^j^sp) regions for three benchmark K values, K = 
0.8, 0.4, 0.1, corresponding respectively to NLSP=fi, bi, Vr- Points below the lines are cosmo- 
phenomenologically allowed. For K = 0.8 we find that, for 32 < tan/3 < 38, the coannihilation 
region is contiguous to the A pole direct annihilation region, producing high values for Ajvlsp 
which are however unrelated to the efficiency of the coannihilation processes. 

parameter K. We will follow here the formalism of ref. [HHI in order to study the loga- 
rithmic dependence of the neutralino relic density . There, the amount of fine-tuning was 
described through the overall logarithmic sensitivity 




where Oj are generic SUSY input parameters. As outlined in JH]; the coannihilation re- 
gions tend to have higher fine-tuning due to the steep rise of the cross sections at the 
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rriQ = 2750 GeV 
mo = 4400 GeV 




Figure 18: The fine-tuning parameter A^{K) = for the four benchmark values 

of mo = 550, 1760, 2750, 4400/GeV. The outher SUSY input parameters are fixed at 
tan/3 = 38.0, M1/2 = 1100/GeV, Aq = 0, hence the benchmark mo correspond respectively 
to (mo/Mi/2) = 0-5, 1.6, 2.5, 4 (see figW^. At the left end of the lines corresponding to 
mo = 2750, 4400 the respective NLSP becomes tachyonic, as indicated. 



onset of coannihilation processes. In the present case we concentrate only in the contri- 
bution A^{K) coming from the parameter K, and study four benchmark cases at fixed 
tan/3 = 38.0, M1/2 = 1100 GeV and Aq = 0. Motivated by figH we choose to plot the 
following values for (mo/Mi/2): 0.5, corresponding to the ordinary stau NLSP coannihila- 
tion region; 1.6, intersecting the low K sneutrino coannihilation region, where fine-tuning 
of K is expected to be reasonably low; 2.5 and 4.0, values which cut the steep part of 
the new coannihilation branch, corresponding respectively to tau sneutrino and sbottom 
NLSP, and where fine-tuning is expected to be rather high. FigUHl reproduces our results: 
for (mo/Mi/2) = 0.5 and 1.6 we get, as expexted, low values for A^{K) < 3. On the 
other hand, the sharp onset of coannihilations for (mo/Mi/2) = 2.5,4.0 entails dramatic 
K fine-tuning peaks, corresponding to the values of K such that m-j^ ~ m/j^ and ~ mi_^ 
respectively. In the case of the sbottom, as suggested by figlH the steepness of (K) 
renders the fine-tuning needed to get ~ particularly large. 

To conclude, we find that the extended sfermion coannihilation channels of mNUSM mod- 
els are characterized by a large amount of fine-tuning in the non-universality parameter K 
in the narrow upper part of the coannihilation branch depicted in figEl Low fine-tuning of 
K is instead required in the lower part of the branch, a.t K < 0.2 and 1 < (mo/Mi/2) ^ 2 
and in the standard stau coannihilation strip. Nonetheless, we stress that these scenarios, 
though motivated in the context of SU (5) SUSY GUTs, are ad hoc sketchily simplified. 
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reducing to an only parameter K the scalar non- universality variables, in order to fo- 
cus on peculiar coanniliilation phenomena. In "realistic" models one could expect to 
reproduce more naturally, hence with a smaller amount of fine-tuning, the new outlined 
coannihilation regions. 

7 Conclusions 

In this paper we studied the cosmological and phenomenological consequences of top- 
down b-T Yukawa coupling unification with minimal non-universal boundary conditions 
in the sfermion masses, inspired by SU (5) GUT. We showed that the /i > case is ruled 
out by the 6-quark mass bound. We stress that this result holds also in the particular 
case of full universality (CMSSM). As for the n < case, b-T YU is compatible with 
the set of all known cosmo-phenomenological constraints, among which the recent results 
from WMAP, for 31 < tan/3 < 45. A large parameter space region is also found to be 
consistent with the 2-a range of determined from r decay data. Further, if one resorts 
to a super-conservative approach [US] to (5a^, the whole allowed regions discussed would 
fulfill the resulting bound. 

We found that the SUSY spectrum allows for new types of coannihilations, namely 
neutralino-sbottom and neutralino-tau sneutrino-stau, that we analysed in detail. We 
fixed three benchmark scenarios for the three possible coannihilation patterns, including 
the CMSSM-like case of neutralino-stau, and we showed for these three cases the behavior 
of the cosmological and phenomenological constraints as functions of the LSP mass. We 
then discussed the cosmologically allowed regions for the two types of new coannihilations 
at various values of tan (3, and the main channels contributing to the neutralino relic den- 
sity suppression. An analytical approximate treatment of these channels is given in the 
Appendix. 
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A Neutralino relic density calculation in the pres- 
ence of coannihilations 



The starting point to compute the neutrahno rehc density is the generahzation of the 
Boltzmann equation to a set of coannihilating species [101 ^] 

^ = -3Hn - {a,sv) {n^ - nl^) , (32) 

where H is the Hubble constant, n = Xlili is the total number density summed over 
all coannihilating particles, n^q is the equilibrium number density, which in the Maxwell- 
Boltzmann approximation, valid in the present cases, reads 



rp N 



(33) 



Qi being the internal degrees of freedom of particle i of mass m,, T the photon temperature 
and K2{x) a modified Bessel function. In eq. ()32|1 aefr is the effective cross section, defined 

as 

N 

dcff = ^ aij Ti Tj. (34) 

In its turn, aij is the total cross section for the processes involving ij (co-) annihilations, 
averaged over initial spin and particle- antiparticle states. The coefficients r,, in the rea- 
sonable approximation where the ratio of the number density of species i to the total 
number density maintains its equilibrium before, during and after freeze out, are defined 
as 



<, ^ Qi (1 + A,)^^/^-^'- 

'^eq 5'tot 



(35) 



where 



TV 

= Y.9.il + ^.f" e-^^ ^ A. = ^^^i^, X ^ ^. (36) 
1=1 

From ()35j) it is apparent that only species which are quasi degenerate in mass with the 
LSP can effectively contribute to the coannihilation processes, since large mass differences 
are exponentially suppressed. Once numerically determined the freeze-out temperature 
Tp ^2); the LSP relic density at the present cosmic time can be evaluated by ^T] 

gj Mp Xp (Teff 

where Mp = 1.22 x 10^^ GeV is the Planck scale, ~ 81 is number of effective degrees 
of freedom at freeze out, xp = m^/Tp and 

/■oo 

aeff = Xp (o"efrf . (38) 

J Xp 
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keeps track of the efficiency of post-freeze-out annihilations. 

In many cases, one can approximate the thermally averaged product of the relative velocity 
and the cross section of the (co-) annihilating particles through a Taylor expansion 



This approximation is not accurate near s-channel poles and final-state thresholds, as 
pointed out in ref. fn\ I13j. In all other cases, and namely in the great part of the ones 
studied in the present paper (an exception is given in figE)), one can proceed and calculate 



where the sum is extended to all the possible pairs of initial sparticle states, and the 
coefficients 



We list in Tables El -El all the possible annihilation and coannihilation processes involving 
the neutralino, the stau, the sbottom and the tau sneutrino. Only a subset of these 
reactions effectively contribute to the reduction of the neutralino relic density, as described 
in sec J6.1l and l6.2l In particular, contrary to the case of neutralino annihilation, the largest 
contributions to ()40|) arising from coannihilation processes come from the aij coefficients. 
In sec JA.2l we give the analytical form of the aij for the "new" coannihilations which arise 
in the present context of mNUSM. Namely, we study the most relevant processes in the 
cases of sbottom-neutralino, sbottom-sbottom and stau-tau sneutrino coannihilations. 
A numerical check of the formula given in sec JA.2l consisting in a comparison between 
the computation outlined in this Appendix and the numerical result given by micrOMEGAs 
confirmed the expected validity of this approximate treatment to a satisfactory extent, in 
the regimes not affected by the A-pole effects. 



(39) 




(40) 




(41) 
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A.l The coannihilation processes 

We report in Tables El - El the complete list of all (Co-)annihilation processes involving 
the lightest neutralino, sbottom and stau as well as the tau sneutrino. For a given couple 
of initial sparticles we list both all the possible final states and the tree level channels 
relative to any final state, c means ^-pO'i^Hcles contact interaction, while s(X), t{X) and 
u{X) mean an s, t or u channel where X is the exchanged (s)particle. d and u indicate 
respectively the down and up-type quarks, while I and z/ the charged leptons and the 
neutrinos of any family, where not differently specified. / stands for a generic fermion 
(quark or lepton) 



Initial state 


Final states 


Tree level channels 


X X 


hh, hH, HH,AA, ZZ,AZ, 


s{h), s{H), t(x°), u(x°) 




h[H] A, h[H] Z 


s{A), s{Z), t(x°), u(x°) 




W+ H+ H- 


s{h), s{H), s{Z),tixj), n{xj) 






s{h), s{H), s(A),t(x7), u(x7) 




fl 


/ 1 \ /tt\ /4\ / r- 7 \ , / P \ / l>\ 

s(/i), s(/7), s(A), s(Z), t(/), u(/) 


X bi 


b h[H], b Z 


s(&), t(6i,2), u(x°) 




bA 


s(6),t(62), u(x°) 




bl,b g 


s(6), t(6i) 




t t W- 


S(&), t(ti,2), U(X7) 


X n 


T h[H], T Z 


s(r), t(fi,2), u(x°) 




tA 


s(r),t(f2), u(x°) 




T 7 


s(r), t(fi) 






s(r), t(z>^), u(x7) 


X 


Vr h[H], Vr Z 


S(Z/^), t{Ur), U(X-) 




Vr A 






T H+, T W+ 


S(z^r), t(fi,2), U(X7) 



Table 2: Neutralino annihilations and coannihilation with sbottom, stau and sneutrino. 
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Initial state 


Final states 


Tree level channels 


h h 


h h 




k k 


hh, hH, H H ,Z Z 


s{h), s{H), t(6i,2), u(6i,2), c 




A A 


s(h), s(H), tih), U(b2), c 




Z A 


s{h), s{H), t{h), u{h) 




A h[H] 


s{Z), tih), nih) 




Z h[H] 


S{Z), t(6i,2), U(6i,2) 




W+ W-, H+ H- 


s{h), s{H), s{Z), t{h:2). c, s(7) 






s{h), s{H), t{h,2) 




tt,uu 


s(h), s(H), s(Z),s(7),s(y),t(x7) 




hh 


s{h), s{H), s(Z),s(7),s(^),t(x7),t(^) 




dd 


s{h), siH), s{Z), s(7), s(^) 




11 


s(/i), siH), s(Z),s(7) 




V V 


/ r7\ 

s{Z) 




7 7, 7 Z, 7 g, Zq 


t(6i), u(6i), c 




^h[Hlgh[H] 


t(6i), u(6i) 




9 9 


s(^), t(6i), u(6i) , c 




hr 


t(x°) 




h T 






t Vt 






h Vt 


t(x?) 




t T 


t(x7) 




h Vj- 


t(x°) 



Table 3: Sbottom annihilations and coannihilations with stau and sneutrino. 



30 



Initial state 


Final states 


Tree level channels 




T T 


t(x?), u(x°) 




/i /i, /i //, H H ,Z Z 


s(/i), S{H), t(fi,2), U(fi,2), C 




A A 


s(/l), S(//), t(f2), U(f2), C 




A Z 


s(/l), S(//), t(f2), U(f2) 




A h[H\ 


S(Z), t(f2), U(f2) 




Z h[H] 


S(Z), t(fi,2), U(fi,2) 




W+ H+ H- 


s(/i), s(if), s(Z), u(z>^), c, s(7) 






s(/i), s(iJ), u(z>.) 




u u, d d, I I 


s(/^,), s(Z), s(7) 




T T 


s(/i), s(i/), s(Z), t(x°), s(7) 






<Z), t(x7) 




V V 


s(Z) 






C, t(Ti), u(Ti) 




1 "-[-"J 






T 


m, u(x7) 




u djU I 


s{H-), s{W-) 






s{W-), s(i/-),t(x°) 




z w- 


s(Vr-), u(fi,2), t(z>^), c 




7 W- 


s(Vr-), u(fi), c 




Z H- 


s{H'), U(fi,2) . t(z>^) 




7 H- 


s{H-), u(fi) 




W- h, W- H 


S(//-), S(iy-), U(fi,2) , t{i>r) 




W- A 


s{H-), U(f2) 




H- h, H- H 


s(//-), s{W-), U(fi,2) , t(i>0, c 




H- A 


S{W-), U(f2), C 



Table 4: Stau annihilations and coannihilations with sneutrino. 
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Initial state 


Final states 


Tree level channels 










hh, hH, H H , Z Z 


s(/i), s(iJ), t(z/^), u(z>^), c 




A A 


s(/i), s(i/), c 




Z A 


s(/i), s{H) 




Z h[H] 


s(Z), t(i>^), u(z>^) 




W+ H+ H 


s(/l), S(i/), S(Z), t(fi,2), C 






S(/1), S(i/), t(fi,2) 




u u, d d, I I, 


s(/i), s(i/), s(Z) 




T T 


s(/i), s(if), s(Z), t(x°) 






s(^), t(x°) 




A h[Hl vV 


s(Z) 



Table 5: Sneutrino annihilations. 



A. 2 Relevant approximate formulae for the relic density calcu- 
lation in the coannihilation regions 

The a coefficients for the pair annihilation of neutralinos can be readily and completely 
derived from ref |H5 , while those concerning the stau annihilations and coannihilations can 
be calculated from the formula of ref. [82j and fSlll2]- Some instances of a coefficients 
computations can be found in ref. ^Zj, with the corrections given in [70 . We follow here 
the notation set in :17^. As regards the tau sneutrino (Co-) annihilations, the relative a 
coefficients can be readily derived from those of the stau via suitable replacements in the 
couplings (e.g. gfj^f^^h Qv^urh etc.), in the masses of the (s-)particles (e.g. mf^ f^u^) 
and setting the mixing angle between the stau mass eigenstates to zero. We do not give any 
explicit formula for the sbottom-stau and the sbottom-sneutrino coannihilations since even 
in the transition regions where ~ mj>^ ~ m^^ these processes do not give any sizeable 
contribution to the neutralino relic density suppression (namely, their total contribution 
is always less than 0.5%). Therefore we are left with the "new" (co-)annihilation processes 
involving respectively neutralino- shottom, sbottom-sbottom and stau-tau sneutrino. 
We give here the explicit formula for the dominant processes, as discussed in sec Jfj.ll and 
16.21 In the kinematic part of the a-coefficients we neglect the mass of the final Standard 
Model Particles up to the 6-quark included. Since the mass of the final SM particles msM 
appears there as ''t^'sm/^'susyi corrections are in fact always negligible. On the other 
hand, we keep trace of both mi, and m,- if the couplings or the whole amplitudes are 
proportional to them, as it is the case in some of the considered processes. Moreover, in 
the processes involving the sbottom we followed the approximations of ref. JH] neglecting 
the terms in and aw Bar over a mass means that the mass is divided by the sum 
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of the masses of the incoming sparticles. 5* stands for the three neutral physical Higgs h, 
H and A, and ms for the respective masses, Sx = sin6'x, Cx = cos9x and for the other 
symbols we follow the notation of ref. jHI- 

A. 2.1 Couplings 



(J oyiiiuol 


Expression 


„L 




( (at tan^H- ^^ nibNi^ ^ 


) 




V V 6 / mw cos p y 




-92 ( 
V2 [ 


rribNi^ 2 
Cfc + s^tan^^y 
mvi/ cos p 6 


A^ii j sign(m^ 




9bibig~9bbg=9ggg 


-9s 


^bibg 


V2gs si 


^bibg 


-V2gs Cl 


0® 
bibXi 


2 {\bx. +^£bx>, 




^blbx^ 


2 {^hbx, ~ ^iibxi, 




9bibigg 


9l 


9h[H]fifi 


92 

2cw 
1 92mr 

' MwCi3 


mzSa+f3[-Ca+p\ ((1 " 2s^)4 + 2s^c?^ 

'^rSa[~Co] + Cf ^A^Sq, [— Cq,] + /iCQ-[s 


+ 

.])] 


Q A ■T-i 'Ti 





9WUrfl 


92Sr 

v/2 


9WUrf2 


92Cf 

v/2 


9H+h[H]W- 


92 r 1 

- — Ca-p[Sa-p\ 


9H+AW- 


92 

2 


9H+Vrri 


g2MwSf 
v/2 



Table 6: Relevant couplings used in the formulcE of sec \A.2.^ (I) 
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g Symbol 


Expression 


9h[H]W+W- 


92MwSf3-a[Cf3-a\ 


Qaw+w- 





9h[H]fiT2 






Ar tan/3 +u 

92 "^T 

^-i-y-L \/\/ 


i/Vy ZjUtT 1 


92^W „ 


9wwz 


fl'2 Cw 


9wtb 


92 

V2 


9H+tb 


K cot/?) 


9H+tb 


j'^ — (rrif, tan/3) 


9^ -0 


V2Mvt/ C/3 \ ch/ / 


9- -0 




-0 





R 




L 


92 Uj2m^ 

Cf 

V2Mw Cfs 




g2UjiSf 




92 U*2mr 
V2Mw cp 


^!^rrx, 


-92Vji 



Table 7: Relevant couplings used in the formul(B of sec VK.2.^ (II) 
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A. 2. 2 X — &i coannihilations 




A. 2. 3 6i — 6i annihilations 



54 71 ^m? + ^^^^ 




+ 



2 TT ( m| + m~ ) ( ml + m- ] 



4 ^x. ^x. + (^g^. J (^g^. j' 



^939 ^fifeig + ^bl6i99 (^^bifeigs ^'ftifeig) 



1728 TT m? 
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A. 2. 4 fi — Ut coannihilations 



128 TT mp^mfi 



2 5'Sflfl S'tyt^rTl 



2 Qh+sw- gH+u-rfirrii,^ Qwurfi 9s w+w- {{m^r 



1 - m|± ' M-^ (1 - M2^) 

2 gSfif2 gWyrT2 ^i>r 



1 /2 

([1 - {Mw + Mz)'] [1 - [Mw - Mzf]) 



X 



32 TT mi>^ 



fl'W'Zi>r-fl + Qwvrfi 9wwz ( 1 



+ M| - 1 
2MzMw 



+ 4 g^H^ziJ^fi gworfi gwwz 1 



M 



.f2 \ 2 



16 TT my^nif^ 



2 2 ^bf^t 

QWvrf-i 9wtb ~Q4. 



+ 



9H+tb\'^ + \9H+tb\'^ ■>-nbmt+ (y^+tfe &f +t6 ) ("^? - l) ) 



9H+Vrfl 9Wi>rfl 9Wtb 



[pH+tb^t (1 - "^?) - +t6"^6 (1 + ml) 
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32 TT mp^mfj 



4 

«,fe=i 



m- Re 



4 R-e 



+ E 

i,fe=l 



of .0 fo^' -o) -0 (g^ ~o) 



2 m^^ Re 



^r-iVrXj 



9- -) 9^ - (9^ -) 



2 m?^ Rc 



E 



9~ 



4 2 mf.mj^^Re f fof of -0 9? -0 {9?" —1 1 

i=i j=i 



mi- + rrii^m 

Xi 



Tl 



i=i j=i 



O? 

^riTXi 
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